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Azoreductases are well known for azo pro-drug activation by gut flora. We show that azoreductases have a 
wider role in drug metabolism than previously thought as they can also reduce and hence activate 
nitrofurazone. Nitrofurazone, a nitroaromatic drug, is a broad spectrum antibiotic which has until now 
been considered as activated in bacteria by nitroreductases. The structure of the azoreductase with 
nitrofurazone bound was solved at 2.08 A and shows nitrofurazone in an active conformation. Based on the 
structural information, the kinetics and stoichiometry of nitrofurazone reduction by azoreductase from 
P. aeruginosa, we propose a mechanism of activation which accounts for the ability of azoreductases to 
reduce both azo and nitroaromatic drugs. This mode of activation can explain the cytotoxic side-effects of 
nitrofurazone through human azoreductase homologues. 

Azoreductases are a group of NAD(P)H dependent flavoenzymes that have been identified in a range of 
bacterial species found in the gut including E. coli 1 and P. aeruginosa. 2 Members of this family are also 
found in humans where they are known as NAD(P)H quinone oxidoreductases. 3 The importance of 
azoreductases in drug metabolism in the gut is well known for azo drugs used to treat inflammatory bowel disease 
e.g. balsalazide. 4 Data have shown that these enzymes can reduce a number of other classes of substrate including, 
quinones 1, 5 and nitroaromatics. 6 As a result of the diverse substrate specificity profile it is suggested that the role 
of azoreductases in drug metabolism may be wider than originally thought. 

Nitroaromatic groups are found in a range of drugs currently on the market including the anti-inflammatory 
nimuselide and the parkinson's disease treatment tolcapone. A common feature of some nitroaromatic drugs is 
hepato toxicity resulting from reduction of the nitro group. 7 Nitrofurazone was selected as an example of a 
nitroaromatic drug to study the mechanism of nitroreduction by azoreductases and their human homologues 
NAD(P)H quinone oxidoreductases. Nitrofurazone (Fig. 1) is a topical antibiotic mainly used to treat burns and 
skin grafts. Nitrofurazone is part of the 5-nitrofuran family, which are potent antibiotics, active against both 
Gram-positive and Gram-negative bacteria. Members of the 5-nitrofuran family are being developed as therapies 
for multidrug resistant bacteria e.g. M. tuberculosis. 8 Nitrofurazone and the other members of this family have a 
common mode of action that requires intracellular activation by nitroreductases. As a result, mutations within 
two nitroreductases known as nitrofuran sensitivity genes A and B (NfsA and B) from E. coli provide nitrofur- 
azone resistance. 9 As humans lack nitroreductases such as these, the majority of nitroreduction will be carried out 
by NAD(P)H quinone oxidoreductases. 10 Activation by nitroreductases results in four electron reduction of the 
nitro group to a reactive hydroxylamine metabolite 11 that in combination with the nitroso oxidation product 
causes DNA damage. 12 

Three azoreductases from P. aeruginosa, paAzoRl, paAzoR2 and paAzoR3, have been extensively charac- 
terised. 2 ' 5 ' 13 ' 14 A recent structural study has provided insight into the mechanism by which azoreductases activate 
azo drugs, and also explains their ability to reduce both azo and quinone substrates. 13 The aforementioned ability 
of azoreductases to act as nitroreductases indicates they may also reduce nitrofurazone. 6 ' 15 

The structure of E. coli NfsB bound to nitrofurazone has been solved 11 the structure showed nitrofurazone 
bound in what the authors described as an inactive conformation within the enzyme active site. In order for their 
proposed mechanism of reduction to be correct nitrofurazone has to reverse its binding orientation as a result of 
FMN being in the reduced state. The work described in this manuscript demonstrates that azoreductases can also 
reduce nitrofurazone and describes the kinetics of nitrofurazone reduction by paAzoRl as well as the structure of 
nitrofurazone bound to paAzoRl in an active conformation. Based on this structure, a putative mechanism is 
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Figure 1 | Structure of nitrofurazone. Figure showing the two nitrofurazone tautomers that occur in equilibrium in solution. 



proposed for the reduction of nitrofurazone by azoreductases that is 
also applicable to the mechanism of reduction of nitrofurazone and 
other nitroaromatic drugs by bacterial azoreductases and by their 
human homologues. 

Results 

Nitrofurazone is reduced by paAzoRl. Nitrofurazone was shown to 
be reduced by paAzoRl, apparent K m and V max values were obtained 
via Lineweaver-Burk plots. The apparent K m for paAzoRl was 15.7 ± 
1.6 |iM and the V max was 9.8 ± 0.6 |iM.s _1 .mg of enzyme" 1 . Via 
equations (1) and (2) the rates of nitrofurazone reduction 
and NADPH oxidation were obtained simultaneously under the 
conditions described in the materials and methods. The 
simultaneous rates obtained during nitrofurazone reduction were 
3.12 ± 0.06 |iM.s _1 .mg of enzyme -1 for nitrofurazone and 5.76 ± 
0.25 |iM.s _1 .mg of enzyme -1 for NADPH this gives an estimate for 
the molar ratio of the reaction of approximately 2 moles NADPH 
oxidation per mole nitrofurazone reduction. Mass spectroscopy data 
obtained from a completed reaction carried out in ammonium 
bicarbonate gave a major peak at 184 Da (Supp. Fig. 1). 

Structure of paAzoRl binding nitrofurazone. The structure of 
paAzoRl binding to nitrofurazone was solved at a resolution of 
2.08 A with final R free and R work values of 20.4% and 16.1% 
respectively and with good stereochemistry (Table 1). No 
significant conformational changes were observed between the 
structures of paAzoRl bound to methyl red and paAzoRl bound 
to nitrofurazone (RMSD on carbon ot positions 0.23 A). No 
significant movements of the amino acid side chains within the 
active site were observed when the structures were compared. 

Clear positive difference density was visible for the nitrofurazone 
bound within the active site after a single round of refinement in 
PHENIX (Fig. 2A). Nitrofurazone was built into the density initially 
in one orientation, however upon further refinement it became clear 
there was density for nitrofurazone in two overlapping orientations. 



In both orientations the nitro group of nitrofurazone is positioned 
above the N5 of FMN in contrast to the orientation observed in the 
structure of NfsB bound to nitrofurazone where the nitro is 7.7 A 
from the N5 (Fig. 3F 11 ). Nitrofurazone binds in an overlapping con- 
formation to that previously observed for balsalazide bound to 
paAzoRl (Fig. 2B 13 ) and methyl red bound to the Y131F mutant of 
paAzoRl, 14 as a result it makes many of the same interactions. The 
furan rings of both conformers overlap the salicylate ring of balsa- 
lazide in the corresponding structure and the polar tail of nitrofur- 
azone overlaps with the position of the azo bond of balsalazide as well 
as partially overlapping the 4-aminobenzoyl-P-alanine ring of balsa- 
lazide (Fig. 2B). As the furan ring forms weaker n — n stacking 
interactions, than the benzene rings of balsalazide, the furan ring 
binds at an angle to the plane of the FMN. As well as forming n-n 
stacking interactions with FMN, nitrofurazone forms hydrophobic 
interactions with Tyrl31, PhelOO, Phel51 and Phel73, while the 
nitro group points into a water filled channel close to the side chain 
of Asn99. 

Discussion 

paAzoRl binds nitrofurazone with a significantly lower K m than azo 
compounds 14 while the rate of reduction is also lower. The slow rate 
of reduction is likely to be due to competition between nitrofurazone 
and NADPH for binding to the active site. Due to the very high iC m of 
NADPH (—1.2 mM 14 ) the tight binding of nitrofurazone causes 
substrate inhibition. Apparent V max is higher than that previously 
observed in R. sphaeroides azoreductases (0.23 |iM.s _1 .mg of 
enzyme" 1 ). 16 The rate of nitroreduction by paAzoRl is significantly 
slower than in nitroreductases such as NfsB (/c cat 2.3 s _1 for paAzoRl 
compared to 225 s" 1 for NfsB 11 ) suggesting that in bacteria the 
majority of activation is carried out by nitroreductases. 

The calculated molar ratio of approximately 2:1 for the rate of 
NADPH oxidation to nitrofurazone reduction is in-line with the data 
that were obtained from Electrospray ionisation, time of flight mass 
spectroscopy (ESI-TOF-MS Sup. Fig. IB). The mass spectroscopy 
data agree with previously published data 11 that show the final reac- 
tion product is consistent with the hydroxylamine reductive meta- 
bolite of nitrofurazone (MW =184 Da). 

As when determining the mechanism of reduction of azo drugs, 13 
the tautomeric state of nitrofurazone bound to the enzyme is import- 
ant, the tautomer that binds paAzoRl can be determined from the 
structure. Nitrofurazone can exist in two tautomers (Fig. 1) that are 
broadly equivalent to the azo and hydrazo tautomers of balsalazide. 
Based on previous work on this enzyme, one might expect that the 
azo form would be the more likely substrate as it has a more quinone- 
like ring. 13 The azo form of the compound would be planar due to its 
extensive conjugated system, the density for nitrofurazone in the 
structure however is not planar (Fig. 2A) and this is more in keeping 
with the hydrazo tautomer where one of the nitrogen atoms is sp 3 
hybridised. This hypothesis is in agreement with the proposed con- 
formation of the bound balsalazide with where the "azo" bond is in a 
similar conformation (Fig. 2B). 

In the structure of NfsB binding to nitrofurazone, 11 nitrofurazone 
is observed to bind within the active site in a non-planar conforma- 
tion (similar to that observed in the structure with paAzoRl) with its 
nitro group bound 7.7 A from the N5 of FMN. In this structure 
nitrofurazone is anchored within the active site via a network of nine 



Table 1 | Processing and refinement statistics for pa AzoRl binding 
nitrofurazone. 

paAzoRl nitrofurazone 

Space Group P3 q 2 1 

oc, p, y (°) 90, 90, 1 20 

a, b, c (A) 82.20,82.20, 108.46 
Processing Statistics 

Resolution Range (A) 29.76-2.08 
Unique reflections 0 25,761 (1 ,794) 

Emerge 0 0.061 (0.766) 

<l/a(l)>° 18.6(2.0) 

Completeness %° 99.6 (95.3) 

Multiplicity 0 5.4 (5.2) 
Refinement Statistics 

R W ork% 16.07 

Rfree% 20.35 

RMS bond angle (°) 0.98 
RMS bond length (A) 0.007 

a Numbers in parentheses are for the highest resolution shell. 
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Figure 2 | Binding of nitrofurazone to paAzoRl. A) The binding of nitrofurazone in its two alternate conformers are shown with brown carbon atoms. 
The FMN is shown with yellow carbon atoms, both conformations of nitrofurazone are shown with brown carbon atoms. Important side chains are 
labelled. Blue mesh is the refined 2F 0 -F C map contoured at 1 a, while black mesh is unbiased F 0 -F c map at 2.2a. The black sphere represents a bound water 
molecule. (B) Comparison of balsalazide and nitrofurazone binding to paAzoRl. Colouring of nitrofurazone is as in (A) and balsalazide is in purple. This 
figure was produced in Pymol vl.l. 31 



short range hydrogen bonds to a combination of backbone amides, 
the ribose chain of FMN and the side chains of Glul65 and Asn71. 
This is in sharp contrast to the entirely hydrophobic interactions 
nitrofurazone makes within the active site of paAzoRl. As the bind- 
ing conformation observed in the structure with NfsB is incom- 
patible with the proposed mechanism of reduction it was 
speculated that when binding to the reduced enzyme, nitrofurazone 
binds in the reverse orientation, which would mean it is in a similar 
orientation to that observed in the paAzoRl -nitrofurazone strucu- 
ture. High resolution structural data (1.9 A) from a closely related 



nitroreductase (88% identical - all active site residues are conserved) 
from E. cloacae in both the reduced and oxidised states show that 
reduction of FMN has little effect on surrounding active site resi- 
dues. 17 Lower resolution (2.5 A) structures of the reduced and oxi- 
dised forms of NfsB are also available and show there is no significant 
conformational changes in the active site residues as a result of FMN 
reduction. 18 It is therefore unlikely that reduced NfsB would bind 
nitrofurazone in a different conformation. The large number of 
interactions which are involved with the positioning of nitrofurazone 
in the conformation observed in the oxidised NfsB structure would 




Figure 3 | Comparison of the overall structures (A-C) and nitroaromatic binding (D-F) of paAzoRl (A,D), hNQ02 (B,E) and NfsB (C,F). (A-C) The 
overall structure of the dimeric forms of each of the three proteins. In all three panels monomer A is in green and B is in turquoise. (D-F) A detailed view 
of nitroaromatic drug binding within the active site of the enzymes (nitrofurazone in paAzoRl and NfsB and CB1954 in hNQ02). In all panels the 
flavin cofactor is in yellow (FMN for paAzoRl and NfsB, FAD for hNQ02). The Nitroaromatic drug is in brown in panels D-F. The structure of hNQ02 
bound to CB1954 is from PDB 1XI2 20 and the structure of NfsB bound to nitrofurazone is from PDB 1YKI. 11 
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mean that most likely alternate conformations will be energetically 
less favourable. 

The proposed mechanism for nitroreduction, based on the azo 
tautomer being reduced, is shown in Fig. 4. Reduction reactions 
carried out by azoreductases proceed via a Bi-Bi ping-pong mech- 
anism 1, 14 therefore as occurs during azoreduction, 13 substrate bind- 
ing is preceded by NADPH binding. While bound to paAzoRl, 
NADPH transfers a hydride ion to the N5 of FMN which becomes 
FMNH". NADPH is subsequently displaced by nitrofurazone which 
accepts the hydride from FMNH" (Fig. 4A). The position of the 
nitrogen of the nitro group of each conformer of nitrofurazone 
within 3.6 A of the N5 atom of FMN suggests that this is the target 
for hydride transfer. As the nitrogen atom of a nitro group is posi- 
tively charged this also makes it a good candidate for transfer of the 
negatively charged hydride. To complete the reduction reaction a 
proton is accepted from the water molecule (Fig. 4 A) which is stably 
bound to the protein within 3 A of the nitro group of nitrofurazone 
(Fig. 2A). The reduction takes place twice according to the same 
mechanism, in order to get complete formation of the hydroxyl- 
amine species, therefore, a double Bi-Bi ping pong mechanism can 
also be proposed for nitrofurazone reduction in paAzoRl. 

In addition it was noted that in the structure of paAzoRl binding 
balsalazide, the oxygen of the quinoneimine tautomer of balsalazide 
is —2.7 A from the ketone of the Asn99 side chain. In the structure of 
paAzoRl binding nitrofurazone one of the oxygen atoms of the nitro 
group is a similar distance from the carbonyl of Asn99. This posi- 
tioning would allow hydrogen bond formation after reduction of 
either substrate which may help to stabilise the reduced form of 
the compound and thus may facilitate hydride transfer. 

Nitrofurazone is known to be carcinogenic in mammals. 19 Unlike 
bacteria, humans have no homologues of NfsA and NfsB, as a result 
the majority of nitrofurazone activation in humans is proposed to 
occur via azoreductase-like enzymes. Humans have two azoreduc- 
tase homologues, NAD(P)H quinone oxidorductase (hNQO) 1 and 
2. Although bacterial azoreductases share minimal sequence similar- 
ity with hNQOl and hNQ02 13 they are structurally related (Fig. 3 A 
and B 20 ' 21 ) and carry out many of the same reactions. 3 Both hNQOl 
and hNQ02 have been shown to activate the anti- tumour pro -drug 
CB1954 via nitroreduction. 10 Activation of nitrofurazone by these 
enzymes may result in the carcinogenic properties in mammals 
and may be the cause of contact dermatitis experienced with nitro- 
furazone. The structure of hNQ02 binding CB1954 (Fig. 3E - PDB: 



1XI2 20 ) shows the nitro group that is reduced, binds with the nitrogen 
within 3.6 A of the N5 of FMN in a similar position to the nitro group 
of nitrofurazone bound to paAzoRl. As a result it is proposed that 
other nitroaromatics are reduced by azoreductases via the same 
mechanism as is described here for nitrofurazone. 

Methods 

All reagents were obtained from Sigma- Aldrich (Poole, UK) unless otherwise specified. 

Enzymology. Recombinant paAzoRl was prepared as described previously. 2 Unless 
otherwise stated assays were carried out in 96-well format in a reaction volume of 
100 uL using a Fluostar Omega plate reader from BMG labtech (Aylesbury, UK) 
typically using 10 jo.g paAzoRl per well, 750 uM NADPH and varying concentrations 
of nitrofurazone in 20 mM Tris-HCl pH 8.0, 100 mM NaCl. 

In order to estimate the stoichiometry of the enzymatic reduction of nitrofurazone 
by paAzoRl, the following reaction was carried out: 50 uM nitrofurazone, 400 uM 
NADPH, 10 ug paAzoRl. Both the oxidation of NADPH and the reduction of 
nitrofurazone were simultaneously followed spectrophotometrically at 340 nm and 
405 nm respectively. The molar extinction coefficient of each species was determined 
at each A, max (s N adph, 340 3,300 M 'cm \ s NA dph, 405 0 M 'cm ', SNitrofurazone, 405 
10,800 JVT'cm" 1 and SNitrofurazone, 340 9,500 M _1 cm _1 ). Initial velocities were esti- 
mated via the following equations: 

equation (1) v i(340 nm) ^ 



( d[NADPH] 
X dt " 



ENADPH, 340 + 



dt 

d[nitrofuraz one] 
Ft 



^■nitrofuraz one, 340 



equation (2) v i(405 nm) = 



dOD (405 nm) _ g ( d[nitrofuraz one] 



&nitrofuraz one, 405 



dt V st 

As s nadph 405 is 0, it is straightforward to determine the rate of nitrofurazone 

dlnitrofurazone] , , dlNADPH] 
reduction — . This value is used to estimate — in 

, x _ r dlNADPH] dlnitrofurazone] 
equation (1). The ratio or — to — provides an estimate tor 

the ratio of NADPH oxidation to nitrofurazone reduction. 

Identification of the products of nitrofurazone reduction. Identification of the 
final reduction products was carried out via ESI-TOF-MS. The reduced nitrofurazone 
was produced via the following reaction: 25 uM nitrofurazone, 0.5 mM NADPH and 
100 ug paAzoRl in a final volume of 1 mL containing 25 mM ammonium 
bicarbonate pH 7.9 (<0.5 mM NaCl). The reaction was monitored at OD 405 on a 
Hitachi (Mannheim, Germany) U-2001 spectrophotometer and allowed to 
proceed to >98% completion. paAzoRl was removed via a Sartorius (Ashtead, UK) 
3 kDa molecular weight cutoff vivaspin centrifugal concentrator. The reaction 
mixture containing reduced nitrofurazone was lyophilised with a Lab mode (Milton 
Keynes, UK) Scan Vac, before resuspension in 1 mL methanol. ESI-TOF-MS was 
carried out on a Waters (Elstree, UK) micromass mass spectrometer in negative 
ion mode. 
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Figure 4 | Mechanism for nitrofurazone reduction. (A) First round of nitrofurazone reduction. (B) Second round of nitrofurazone reduction resulting 
in hydroxylamine formation. The proton in both (A) and (B) is likely to be donated by a stably bound water molecule observed in the structure (Fig. 2A). 
Before (A) and between (A) and (B) NADPH binds to paAzoRl and transfers a hydride to FMN before being replaced by nitrofurazone. 
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Structure determination. paAzoRl was mixed with nitrofurazone to give a protein 
solution of 23 mg.mL" 1 paAzoRl in 20 mM Tris-HCl pH 8.0 with 2 mM 
nitrofurazone and incubated for an hour before being washed with 20 mM Tris-HCl 
pH 8.0 to remove DMSO. Crystals grew isomorphously under the same conditions as 
previously described, 2 briefly 2 uL protein complex was mixed with 2 uL 1.6 M 
(NH 4 ) 2 S0 4 0.1 M HEPES pH 7.5 in sitting drop format. Data for paAzoRl- 
nitrofurazone were collected at beamline 102 at the Diamond Light Source (Oxon, 
UK) with a Quantum ADSC detector. Data were integrated via XDS, 22 and scaled and 
merged with SCALA 23 both were run as part of xia2. 24 The structure was solved with 
the program PHASER, 25 using the structure of paAzoRl stripped of its ligands (PDB: 
2V9C 2 ) as the search model. The atomic model was rebuilt and refined with Coot, 26 
Refmac 5.5 27 and PHENIX. 28 Initial translation, libration and screw (TLS) parameters 
were determined with the program TLSMD 29 and TLS refinement was performed 
with PHENIX. Restraints for nitrofurazone were generated in eLBOW. Water 
molecules were added with PHENIX and Coot. The final model was validated via 
MolProbity 30 and all residues were shown to lie within the allowed regions of the 
Ramachandran plot (97.8% in the preferred region and 2.2% in the allowed region). 
Data collection and refinement statistics are shown in table 1. The structure was 
deposited in the protein data bank (PDB ID: 3R6W). 
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